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The dynamics at the critical-point of a general first-order quantum phase transition in a finite 
system is examined, from an algebraic perspective. Suitable Hamiltonians are constructed whose 
spectra exhibit coexistence of states corresponding to two degenerate minima in the energy surface 
separated by an arbitrary barrier. Explicit expressions are derived for wave functions and observables 
at the critical-point. 
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Quantum phase transitions have become a topic of 
great interest in different branches of physics [J, M S] ■ 
These are structural changes occurring at zero temper- 
ature as a function of a coupling constant and driven 
by quantum fluctuations. The phase transitions are said 
to be of second-order when the system changes contin- 
uously from one phase to another, and of first-order 
when there is a discontinuous change between two co- 
existing phases. Advanced experiments have identified 
such structural changes in a variety of mesoscopic sys- 
tems, e.g., nuclei, molecules and atomic clusters, where 
the transitions are between different shapes or geomet- 
ric configurations. An important issue in these systems 
is to understand the modifications at criticality due to 
their finite number of constituents. In the present work 
we study this question in connection with nuclei exem- 
plifying a finite system undergoing a general first-order 
shape-phase transition. 

The role of a finite number of nucleons can be ad- 
dressed in the interacting boson model (IBM) which 
describes low-lying quadrupole collective states in nuclei 
in terms of a system of N monopole (s) and quadrupole 
(d) bosons representing valence nucleon pairs. The model 
is based on a U(6) spectrum generating algebra and its 
three solvable dynamical symmetry limits: U{5), SU{3), 
and 0(6), describe the dynamics of stable nuclear shapes: 
spherical, axially-deformed, and 7-unstable deformed. A 
geometric visualization of the model is obtained by an 
energy surface defined by the expectation value of the 
Hamiltonian in the coherent (intrinsic) state [1, 0] 



\P,T,N) = (iV!)-i/2(6t)^|0) 
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where bl = {l+P'^)-^/'^[Pcosjdl+f3smj{4+dl2)/V2+ 
s"^]. For the general IBM Hamiltonian with one- and two- 
body interactions, the energy surface takes the form 



[a(3^ - b(3^ cos 3j + cl3^ 
(I + /32)2 
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The coefficients EQ,a,b,c involve particular linear com- 
binations of the Hamiltonian's parameters 0, @]. The 
quadrupole shape parameters in the intrinsic state char- 
acterize the associated equilibrium shape. For /? > 
the intrinsic state is deformed and represents, in a varia- 
tional sense Q , a ground band whose rotational members 



are obtained by standard angular momentum projec- 
tion. Particularly relevant for the present work are states 
of good 0(3) symmetry L projected from the prolate- 
deformed intrinsic state |/3, 7 = 0; TV) of Eq. ([!]), 
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\f3;N,L,M) cx Pj^\(3)\ ' T'lmI/?, 7 = 0; TV) 

^ fdx[l + (3'P2ix)fPLix) . (3) 
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Here Pl{x) is a Legendre polynomial with L even (0 < 
L < 2N) and r^'(/?) is a normalization factor. The 
projected states with fixed N and L, involve a mixture of 
components \N,nii,T, L) with quantum numbers related 
to the ;7(6) D U{5) D 0(5) D 0(3) chain. In general 
the L-projected states \/3;N,L) interpolate between the 
U{5) spherical ground state, \s^) = \N,na — t — L — 0), 
at /? = 0, and the prolate-deformed ground band with 
SU{3) character {X,fi) = (2iV,0), at /3 = %/2. 

Shape-phase transitions in nuclei have been stud- 
ied in the geometric framework of a Bohr Hamilto- 
nian for macroscopic quadrupole shapes. Analytic so- 
lutions, based on infinite square-well potentials, called 
E(5) [H and X(5) 0, have been shown to be relevant to 
the dynamics at the critical-point of second- and first- 
order transitions respectively, and empirical examples 
have been presented [l^l- Phase transitions for finite 
N can be studied by an IBM Hamiltonian of the form. 
Hi -\- gH2, involving terms from different dynamical sym- 
metry chains The nature of the phase transition 
is governed by the topology of the corresponding en- 
ergy surface which serves as a Landau's potential. 
IBM Hamiltonians of this type have been studied exten- 
sively [H [H, [H, Q [H 0, [O, [H, p, revealing the 
finite-N attributes of the phase transitions. The crit- 
ical t/(5)-0(6) Hamiltonian was found to exhibit a 7- 
independent flat-bottomed energy surface, characteristic 
of a second-order phase transition. The critical U{5)- 
SU{3) Hamiltonian displays an energy surface with de- 
generate minima separated by an extremely low barrier, 
hence corresponds to a very specific (and non-generic) 
first-order phase transition. It is clearly desirable to ex- 
plore, from an algebraic perspective, the finite-N proper- 
ties of a general first-order critical-point with an arbitrary 
barrier. This is the goal of the present investigation. 
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FIG. 1: The energy surface, Eq. Q, at the critical-point of 
a first-order phase transitio n. The position and height of the 
barrier are /3+ = (-1 -f a/1 + /3g )//3o and h = ,f(/3+) = 
(-1 + VTT^ )V4, with /(O) = /(/3o) = 0, /((») = 1. 



In a first-order phase transition the Landau's potential 
has two coexisting minima which become degenerate at 
the critical-point. Such a structure involving spherical 
and prolate-deformed minima can be accommodated in 
the energy surface of Eq. ^ provided a,b > and = 
Aac at the critical-point. For 7 — the corresponding 
critical energy surface can be transcribed in the form 



E,r^{P^^^) - Eo+cN{N -l)f{l3) 



(4) 



Ecri{P,l = 0), shown in Fig. 1, behaves quadratically 
near the minima at (3 = and P = Pq — 2a/b > 0, and 
approaches a constant for large p. The value of Pq deter- 
mines the position (/? = and height (h) of the barrier 
separating the two minima in a manner given in the cap- 
tion, h and /3+ increase as a function of Pq, and hence 
both low and high barriers can be considered. To con- 
struct a critical Hamiltonian with such an energy surface, 
it is convenient to resolve it into intrinsic and collective 
parts 0,11, 



(5) 



Such decomposition can be done exactly for any IBM 
Hamiltonian, as elaborated in great detail in 0, Si- 
Adapting this procedure to the critical Hamiltonian, the 
intrinsic part (Hint) is defined to have the equilibrium 
condensate \ P = Po,^ = 0; N) , Eq. ([!]), as an exact zero- 
energy eigenstate and to have an energy surface with a 
structure as in Eq. Hint is found to have the form 



Hi, 



(6) 



with Pl^{Po) = /3ost4 + ^ (dtdt)J^2)^ p^^i^p^) = 

l)^P2,-/i(/3o)- Its energy surface coefficients of Eq. ^ 
are a = /i2/^0' ^ = 2/i2/3o! c = /12, Eq = 0, hence satisfy, 
for /i2 > 0, the aforementioned conditions of a first-order 



critical-point. Since Hint is rotational-scalar, it has, 
by construction, the L-projected states |/3 — Pq;N,L) 
of Eq. jS]) as solvable deformed eigenstates with en- 
ergy E = Q. It has also solvable spherical eigenstates: 
\N,nd ^ T = L = 0) = \s^) and | AT, = r = L = 3) 
with energy E = and E = 3/i2 [PiiN - 3) + 5] respec- 
tively. For large N the spectrum of Hint is harmonic, 
involving quadrupole vibrations about the spherical min- 
imum, and both P and 7 vibrations about the deformed 
minimum with frequencies e, and given by d 



e/3 = h2 PlN 



Q{l + Pl)-^ep 



(7) 



For the acceptable range < Pq < 1.4, the 7-band 
is expected to be considerably higher than the /3-band. 
All these features are present in the exact spectrum of 
Hint shown in Fig. 2, which displays a zero-energy de- 
formed (K — 0) ground band, degenerate with a spherical 
{ud — 0) ground state. The remaining states are either 
predominantly spherical, or deformed states arranged in 
several excited K = bands below the 7 band. The co- 
existence of spherical and deformed states is evident in 
the right portion of Fig. 2, which shows the rid decom- 
position of wave functions of selected eigenstates of Hint ■ 
The "deformed" states show a broad distribution typi- 
cal of a deformed rotor structure. The "spherical" states 
show the characteristic dominance of single Ud compo- 
nents that one would expect for a spherical vibrator. 

The collective part (He) is composed of kinetic terms 
which do not affect the shape of the energy surface. It 
can be transcribed in the form 0, Q 



Hr 



C3 Co{3) - 6nd + C5 Co(5) - 4ri, 



-C6 



Ef) 



(8) 



where N — hd + fis, na and hs are the total-boson, d- 
boson and s-boson number operators respectively. Here 
Cg denotes the quadratic Casimir operator of the group 
G as defined in 0]. In general, the collective and intrin- 
sic parts of the critical Hamiltonian ^ do not commute. 
Table I shows the effect of different rotational terms in 
He- For the high-barrier case considered here, {Pq = 1.3, 
h — 0.1), the calculated spectrum resembles a rigid-rotor 
{E ^ aNL{L + 1)) for the ca-term, a rotor with cen- 
trifugal stretching {E ~ aNL{L + 1)- bN[L{L + for 
the C5-term, and a X(5)-like spectrum for the cg-term. 
In all cases the B(E2) values are close to the rigid-rotor 
Alaga values. This behaviour is different from that en- 
countered when the barrier is low, e.g.^ for the critical 
U{b)-SU(i) Hamiltonian (corresponding to Pq — l/2\/2 
and h k, 10""^), where both the spectrum and E2 rates 
are similar to the X(5) predictions (T9| . 

To gain more insight of the underlying structure at 
the critical-point, we focus the discussion to the lowest 
coexisting spherical and deformed states. The need for 
angular momentum projection from mean field solutions 
in regions of shape-coexistence is well known, e.g., in 
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FIG. 2: Left portion: spectrum of Hint, Eq. with /i2 ~ 0.1, /3o = 1-3 and A'^ — 10. Right portion: the number of d bosons 
(ud) probability distribution for selected eigenstates of Hint- 



microscopic calculations employing Hartree-Fock Bogoli- 
ubov states d^]. In the present case, the matrix ele- 
ments of the critical Hamiltonian ^ in the states ([3]), 



Ei''\p) = {P;N,L\H,„\P;N,L) 
be evaluated in closed form 



Ei^\f3) +Eq, can 



4^)(/3) = /i2(/3-/3o)'sf2+c3 L{L + 1)-6D 



C5 



(AT) 

2.L 



2,L 



and then examining the 2x2 potential energy matrix, 

K,,{(3) = (*,| I*,) = % + , (13) 

which reads 

Kii = 0, ki2 = ~CQl3''N{l-rl^)-^'^r^2 , 

4=l,(/3) + 2c6/32iVr?J(l-r?2)-\ (14) 



11 — 

k22 - 



+ C6 iV(iV-l)-(l + /32)2 5(;i) 



(9) 



Here 1) 



(TV) 



with iV = A^(A^ ~ !)• The derived eigenvalues of the 

matrix serve as eigenpotentials, £'^^q(/3), and the cor- 

(±) 



denoted by 5'J . Specifically, 



S. 



(N) 



o{N) o(N- 
■l)-2(iV- 



1) yim 



'-'l.L 



and T>^1 denote the expecta- responding eigenvectors, I^^Jq), are identified with the 

ground (O^j'') and excited (0^) L = states. The de- 
formed states 1/3; iV, L) of Eq. ([3]) with L > are iden- 
tified with excited members of the ground-band (LX^ 
with energies given by E^j^\(3), Eq. E2 matrix ele- 
ments between these states can be evaluated in closed 
form. For example, for the E2 operator, T{E2) — 

necessary matrix elements 
for transitions involving the L = 0^,2^, 0^ states are 
Ti = {(3;N,L' = 2\\T{E2)\\f]-N,L = 0) and T2 = 



'1,L ' '^2,L ' ^2,L 

tion values in the states \(3;N,L) of fid, ns{ns — 1), 
fidind — 1) and hsfid respectively. All these quanti- 
ties are expressed in terms of the expectation value of 



iN-l)S\ 



N 



AN) 
'2,L ■ 



(N) 



2,L " 

iV(iV- 

is determined by the normalization factors of Eq. ([3]) 



The quantity S[ ^ itself 



■ N 



^\(/3)/r(,^)(/3) . 



(10) {(3;N,L' = 2\\T{E2)\\s'^), 



It also satisfies the following recursion relation 
[N - L/2){2N + L+1) 



MN) _ 
^l,L — 



+ 4){N - 1) + 3 + (/32 - 2)(1 + f3^)S[''-^^ 



(11) 



obtained from the fact that P2^{Po) of Eq. ^ annihilates 
the states \P;N,L) for f3o — (3. The mixing between the 
coexisting spherical and deformed L — states, |0i) = 
\s^) and 1^2) = \f3;N,L = 0), can be studied by first 
transforming to an orthonormal basis 



To 



(1 



/3x)ri^Li(/3)] 



.(2 



1/2 



PNl[N\Tf{P)Y/^ 



(15) 



1*2) = (1- 



,2^-l/2 



102) - ri 



2 19^1, 



ri2 = (0i|02) = [A^!r(v^-°)(/3)]-i/2 ^ 



(12) 



with X — v^J^X- The parameter [3 in the indicated 
wave functions, energies, and E2 matrix elements, plays 
the role of an effective deformation whose value is chosen 
at the global minimum of the lowest eigenvalue, e\^2o(.P), 
of the matrix Kij{(3) (fT^l) . As is evident from Table I, 
this procedure leads to accurate finite-N estimates of 
observables at the critical-point. The rotational terms 
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TABLE I: Excitation energies (in units of E{2^) = 1) and B(E2) values [in units of B{E2;2^ ^ Oj'" = 1)] for selected terms 
in the critical Hamiltonian, Eqs. (O,®,®, with (i) C3/h2 = 0.05, (ii) C5//i2 = 0.1, (iii) ce/h2 = 0.05, and iV = 10, Po = 1-3, 
E2 parameter x = —V7/2. The entries in square brackets [. . .] are estimates based on the L-projected states, Eq. ([3}, with 
(i) l3 = 1.327, (ii) /3 — 1.318, (iii) /3 = 1.294, determined by the global minimum of the respective lowest eigenvalue of the 
potential matrix, Eqs. (|13|) -(|14 | ). The rigid-rotor and X(5) ^ parameter-free predictions are shown for comparison. 



cs/hi = 0.05 C5//12 = 0.1 ce/hi = 0.05 rotor X(5) 



E{4+) 3.32 [3.32] 3.28 [3.28] 2.81 [2.87] 3.33 2.91 

E{6t) 6.98 [6.97] 6.74 [6.76] 5.43 [5.63] 7.00 5.45 

E{8t) 11.95 [11.95] 11.23 [11.29] 8.66 [9.04] 12.00 8.51 

E{10t) 18.26 [18.26] 16.58 [16.69] 12.23 [12.83] 18.33 12.07 

E{Ot) 6.31 [6.30] 6.01 [5.93] 4.56 [5.03] 5.67 

B(S2;4+^2+) 1.40 [1.40] 1.40 [1.40] 1.46 [1.45] 1.43 1.58 

B(S2;6j^^4+) 1.48 [1.48] 1.48 [1.48] 1.55 [1.53] 1.57 1.98 

B{E2;8i^6i) 1.45 [1.45] 1.45 [1.45] 1.53 [1.51] 1.65 2.27 

B{E2-10+^8t) 1.37 [1.37] 1.37 [1.37] 1.44 [1.42] 1.69 2.61 

B{E2;0t '^2t) 0.003 [0.003] 0.003 [0.004] 0.24 [0.18] 0.63 



in He ^ affect the value of the effective deformation. 
This highlights the importance of the coupling between 
the order-parameter (/3) and the soft- (rotational) modes 
at the critical-point of a quantum phase transition 2li ]. 
The characteristic spectra, discussed above, of these ro- 
tational terms can now be understood from their contri- 
bution to E'j^\f3), Eq. dH), and Kij{f3), Eqs. (HSll-dlll). 

As observed from Eq. (fTT|) . s[^j^ shows an L{L + 1) be- 
haviour near the SU (3) value /3 = \/2 and, consequently, 
the contribution of the ca-term to E^^\(3) is linear in 
L{L + 1) while that of the C5- and Cg- terms is parabolic. 
The rigid- (non-rigid) rotor-like spectrum of the C3- (C5-) 
terms then follows from the fact that they contribute only 
to the diagonal matrix-element K22{(3) in Eq. ((T4)) . The 
cg-term contributes both to diagonal and non-diagonal 
matrix elements, thus controls the mixing which is es- 
sential for obtaining an X(5)-likc spectrum. 

To summarize, we have considered properties of a first- 
order quantum shape-phase transition in a finite system 
from an algebraic perspective. A suitable critical-point 
Hamiltonian was constructed by resolving it into intrin- 
sic and collective parts. The intrinsic part generates an 
energy surface with two degenerate minima separated by 



an arbitrary barrier, and its spectrum exhibits coexis- 
tence of spherical and deformed states. The collective 
part contains kinetic rotational terms which affect the 
rotational splittings and mixing. Spectral signatures for 
the case of a high-barrier have been shown to differ from 
those of a low-barrier. The dynamics at the critical- 
point can be described by an effective deformation de- 
termined by variation after projection, combined with 
two-level mixing of i = states. Wave functions of a 
particular analytic form can be used to derive estimates 
for energies and quadrupole rates at the critical-point. 
The intrinsic-collective resolution constitute an efficient 
method for studying shape-phase transitions, since the 
derived Hamiltonian is tailored to reproduce a given en- 
ergy surface which, in-turn, governs the nature of the 
phase transition. Although we have treated explicitly the 
IBM with one type of bosons, the tools developed are ap- 
plicable to other finite systems described by similar alge- 
braic models, e.g., the proton-neutron version of the IBM 
for nuclei [3| and the vibron model for molecules [23 |. 
This work was supported by the Israel Science Founda- 
tion. 
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